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Abstract 

Integral expressions are determined for the elastic displacement and stress fields due to 
stationary or moving dislocation loops in three dimensional, not necessarily isotropic, finite 
samples. A line integral representation is found for the stress field, thus satisfying the ex- 
pectation that stresses should depend on the location of the dislocation loop, but not on the 
location of surfaces bounded by such loops that are devoid of physical significance. 

In the stationary case the line integral representation involves a "vector potential" that 
depends on the specific geometry of the sample, through its Green's function: a specific 
combination of derivatives of the elastic stress produced by the Green's function appropriate 
for the sample is divergenceless, so it is the curl of this "vector potential". This "vector 
potential" is explicitely determined for an isotropic half space and for a thin plate. Earlier 
specific results in these geometries are recovered as special cases. 

In the non stationary case a line integral representation can be obtained for the time 
derivative of the stress field. This, combined with the static result, assures a line integral 
representation for the time dependent stress field. 

1 Introduction 

A general formula for the displacement and stress fields generated by a dislocation loop undergoing 
arbitrary motion in an infinite medium was obtained some years ago by Mura (1963). In this 
formulation, the stress field is written as a convolution of the medium's impulse response with 
a source localized along the dislocation loop (i.e. independent of the loop's slip plane). These 
formulae do not apply to the case of finite samples since they are derived assuming homogeneity 
in space, i.e. a dependence of the Green's function in x — x' , a fact that is no longer true in a 
finite sample. The purpose of this paper is to present a generalization of these formulae to finite 
samples, as well as explicit forms for the geometries of a half space and thin plate in an isotropic 
elastic medium. 

Motivation for this work comes partly from an ongoing project that attempts to understand 
recently observed dynamic instabilities for cracks in thin plates (Sharon, Gross and Fineberg, 1996, 
Boudet, Cilibcrto and Steinberg, 1996) in terms of the interaction of the plate's oscillations with the 
crack tip regarded as a continuous distribution of infinitcsimally small dislocations (Lund, 1996). 
An intuitive thought behind this work was that the stress fields produced by stationary or moving 
dislocation loops in a finite sample should not depend on the slip planes chosen, and that indeed the 
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stresses will be continuous through the slip planes, as it is the case in an infinite medium. This is 
rigorously shown here: the expressions found for the stress fields involve only line integrals over the 
dislocation loops. In addition, Gosling and Willis (1994) have emphasized the importance of dealing 
with free surface effects in order to understand the strain relaxation within thin strained layers over 
a substrate such as are employed in the manufacture of high technology devices. They developed 
a line integral representation for the stresses due to an arbitrary dislocation in an isotropic half- 
space, in which the integrand can be constructed by an explicitely given algorithm. The resulting 
integral can be explicitely evaluated for a dislocation half-line. Line integral representations are 
computationaly more efficient than the more easily derived surface integral representations. 

In order to find the explicit form of the line integrals appearing in the stress fields when the loops 
are stationary, one should calculate a "vector potential" , that is, a vector quantity whose curl gives 
a combination of derivatives of the inhomogeneous portion of the Green's function, appropriate 
to the finite sample in question. We have explicitely done this for two examples: an isotropic 
half space and a thin plate. In addition to the calculation of the "vector potential" in these two 
geometries we have explicitely worked out the case of a stationary screw dislocation perpendicular 
to the free surfaces. These two examples can also be obtained independently with the method of 
Eshelby and Stroh (1951), thus providing a check on the algebra. In a non stationary situation 
the displacement velocity field can be written in terms of line integrals, and from this one deduces 
that the time derivative of the stress field can be written in terms of them too. One finally deduces 
that the time dependent stress field can be written in terms of line integrals since one knows that 
the initial condition (stationary case) and its temporal derivative can both be written in terms of 
line integrals. 



2 Displacement fields due to dislocations in finite samples 

Here we derive a formula for the elastic displacement field produced in a finite sample by a moving 
dislocation loop, that is written as an integration over succesive slip planes. 

We consider small displacements U m (x,t) in an homogeneous elastic medium of density p and 
elastic constants Cijki- A general formula for the displacement, velocity and stress fields generated 
by a dislocation loop undergoing arbitrary motion in an infinite medium was obtained some years 
ago by Mura (1963). The displacement field is written as an integration over succesive slip planes, 
where the displacement field has a discontinuity given by the Burgers vector bi . The velocity and 
stress fields are written as a convolution of the medium's impulse response with a source localized 
along the dislocation loop, i.e. these fields are independent of the choice of slip planes. In this 
section Mura's formula for the displacement field is generalized to a finite sample. 

The particle displacement obeys the wave equation: 
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and normal stress-free boundary conditions at the surfaces of the sample: 

C m -^(x> s ,t)n j (x> s ) = 0, (2) 

with x' s a point on the surface and rij(x' s ) its outwardly-pointing normal vector. The medium's 
impulse response is given by a Green's function Gi m (x,x';t — t'), which is the displacement in 
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direction (i) evaluated at (x, t) produced by a localized impulse force in direction (m) applied at 
position x' and time t' . It is the solution of: 

P^POr, x'; t-t')- C ijkl ^^-(x, x'; t-t')= S m S(x - x')6(t - H) . (3) 

We have explicitely exhibited the x and x' dependence because the order in which they apear is 
important. If one defines as a^ m (x, x'; t — t') the elastic stress associated with the Green's function: 

a? m (x, x';t-t') = C m -^—G lm (x, x';t-t'), (4) 

the free surface boundary condition reads: 

af j m {xs,x';t-t')n {xs)=0. (5) 

Note that, since the medium is bounded, there is homogeneity in time but not in space. 

In order to find the displacement field due to a moving dislocation we start, following Mura 
(1963), with the following identity that holds due to a symmetry of the elastic constants tensor 

(Cijkl = Cklij)' 

c mi -^r{x ,t) -^r (x, x ■ t - t ) = C i3 u — (x , t) (x,x;t-t). (6) 

A single dislocation loop is considered to be present. At the position of the loop itself the derivatives 
of the displacement field have singularities, and across the slip plane the displacement field has a 
discontinuity equal to the Burgers vector hi. The identity in Eq. (|^) is integrated over the volume 
of the elastic sample, excluding a thin tube around the dislocation loop, as well as a thin layer 
encompassing the slip plane. Using the reciprocity relation of the impulse response (Poruchikov, 
1993): 

Gi m {x, x'] t — t) = G m i(x' , X] t — t') , (7) 

as well as Eqns. (Q)-(|^), and integrating by parts, the following identity results: 



dS'jCijki {x', t')G m i(x, x 1 ' ;t — t') — J dx' 'Cijki q x i q x i ( x 'i t')Gmi(x, x';t — t') — 

dSfiijuUitf, t')^r(x, x';t-t')-J dx>C tjkl U t (x>, t')^0{x, x';t-t'), (8) 

where the surface integrals are carried out over the tube and layer just introduced, as well as over 
the sample boundary. The first term on the left vanishes on the surfaces of the sample due to the 
free surface boundary condition, Eq. (||), and it also vanishes on the slip plane due to continuity 
of the stress there. The first term on the right also vanishes on the surfaces of the sample due to 
the free surface boundary condition, Eq. (^|), due to reciprocity, but it will give a contribution on 
the slip plane due to the discontinuity of [/, given by the Burgers vector bi. The first terms on the 
left and on the right of Eq. (||) integrated over the small tube (of radius e) around the dislocation 
loop vanish because of the short distance behaviour of particle displacement Ui (Lund, 1988). 
Using the differential equations satisfied by Ui and G m i one gets from Eq. (||): 

^, d 2 Ui(x\t) ;t -t') = bi [ dS' 3 C ijU ^^{x,x';t-t') 

dt 2 J s{v) 3 dx[ 

dx'^P(x, x'- 1 - 0^(^,0 + 6(t - t')U m (x, t') , (9) 
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with S(t') the slip plane at time t'. Integrating over time t' , one obtains the expression for the 
displacement as an integral over the succesive slip planes: 

U m (x,t) = -b l f dt' [ dS'jCijH^^&x'it-t'). (10) 

J-oo Js(t') OX l 

Using reciprocity, this can be written as 

U m (x,t) = -bi[ dt' f dS'ja^ix^x^-t') . (11) 

J-oo Js(t') 

Attention should be paid to the order in which the variables appear in this formula. In an analogous 
way, for the stationary case, the following displacement field is derived: 

U m (x)=-bi f dS' j( r% m (x',x), (12) 

J S' 



with the static Green's function satisfying: 

a jkl —^—G lm (x,x') = —a% m (x,x') = -6 im 5(x - x') , (13) 

OXkOXj OXj J 

and the free surface boundary condition: 

a° m (x s ,x') nj (x s ) = 0. (14) 
Eq. ([ll]) leads to null normal stresses at the free surfaces. Indeed, from Eq. ([!(]) one gets 

(Jpq(x : t) = —biCpqnm I dt I dSjCijklT{ 7- j (x , X ]t £ ) 

J-oa JS(t>) OX n OX l 

f°° f da Gk 
= -biC m dt' dS'i—^&x'it-t'). (15) 
J-00 Js(t>) 3 dx[ 

Since o^ q k (x $ , x' '; t — t')n q (xg) = 0, this implies that <7 pq (xs,t)n q (xs) = 0, as it should. 

3 Stress fields due to dislocation loops in finite samples 
3.1 Stationary case 

Starting from Eq. (fi2|), the stress field produced by the static dislocation loop is: 

a m {x) = -biCpgim / dS'— a% m (x',x) . (16) 

In the infinite medium case, space homogeneity easily leads to a line integral representation (Mura, 
1963). In order to obtain a line integral representation for this stress field in the finite case, we 
add and substract to the integrand the term 

® Gmt i \ 
(*,*), 
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which is easy to turn into a line integral, and we define a tensor: 

f^(x', X ) EE ±-a?™(x>,x) + JL a Gm {x , x) 

which is divergenceless. Indeed, from Eq. (|l3|), 



d , lmf , , , , d . d 



jfl?(x',x) = Sim(g- + gj)6(x -x') = 0. (18) 
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This means that /^(x', x) can be written as the curl of a "vector potential" A l j m (x' , x): 



f^(x',x)=e jrs ^-A^(x',x), (19) 



' dx, 

and it is possible to write the stress field, Eq. (fl6|), as: 

dx 

When x is not in the slip plane S' the following is true: 

— [x,x)-e vjl e vsr / dS s — 



a pq (x) = -6,^ / dS${/#V,aO - — ^(x',x)} . (20) 



dS'— ^(x',x) = e vjl e vsr / dS' s —*% m (x',x) , (21) 



so that we finally have that the stress field can be expressed as a line integral along the dislocation 
loop: 

a pq (x) = biC pqlm f dl' s {e sj ia% m (x',x) - A\ lm (x',x)} , (22) 



as promised, with A\ (x' , x) a "vector potential" to be determined for each specific sample by use 
of Eqns. ( p"7| ) and (|l9|). This means that of^ m , the Green's function for the finite sample, must 
be known. Its determination, however, is of course independent of the particular dislocation loop 
that might be under consideration. 

Eqn. (22) is similar to Eqn. (12) of Gosling and Willis (1994). The difference between those two 
expressions resides in the fact that Gosling and Willis (1994) split the integrand into the Green's 
function for an infinite medium plus "image" forces that annul the tractions at the free surface. 

3.2 Non stationary case 

Since the problem is homogeneous under time translations it is possible to obtain an expression for 
particle velocity, that involves an integral along the dislocation loop only. Indeed, from Eq. (11): 

^{x,t) = bi dt'{^-[[ dS'^ m (x',x;t-t')}- [ dS' j( r? m (x',x;t-t')}. (23) 
9t J-oo of J s{t , } •%(*') 

The first term is zero since the Green's function vanishes at t — t' = ±oo. Since: 

dS'^Cjpgf dl' q V p (x',t'), 
dS(t')/dt> JL(t') 
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where L(t') is the dislocation loop bounding the slip plane S(t'), V p (x', t') is the loop's local velocity 
and ejp q the completely antisymmetric tensor in three dimensions; we then have: 



! ' (7m -(M) = ~bi I dt 1 I dl'^ m (x',x;t- t')e jpq V p (x',t') . (24) 

iL{t>) 



dt 



From this line integral expression for the displacement velocity field, one obtains an expression for 
the partial time derivative of the stress field that involves line integrals over the dislocation loop: 



-El^ln{x,t) = —Cinkmh I dt'j dl'^-a% m {x\ x; t - t')e Jpq V p (x' ,t') . (25) 

Ct J-oo J L(t') OXk 



This result, together with the result for a static loop (considered as an initial condition), mean 
that the time dependent stress field can be written in terms of integration of line integrals over the 
dislocation loop, and it is independent of the choice of slip planes, as it should. Explicitly, we are 
thinking on the generally valid formula: 

ai n (x,t) = a hl (x,t = 0) + J dt'^(x,t') , (26) 

with both terms on the right hand side having a line integral representation. 

A general result that derives from Eq. (25) is that segments of the dislocation loop on the free 
surfaces do not contribute to d<Ti n /dt: in those segments dl 1 x V is parallel to the normal n to the 



free surface, and then the boundary condition satisfied by of™ at the free surfaces cancels that 
contribution 



4 Examples of explicit "Vector potentials" 
for the stationary case 

In two geometries, a half space and a thin plate, and considering an isotropic elastic medium we 
have determined explicitly the form of the "vector potentials" appearing in the general line integral 
representation of Eq. (|22j ) for the stress field generated by a dislocation in a finite sample. The 
derivation will be explained in the following, with further details to be found in Appendices I and 
II. 

The elastic medium fills the half space z > 0, and the thin plate the space between —h < z < h, 
where z is the coordinate whose direction is perpendicular to the free surfaces of these geometries 
(we will call R a coordinate in the planes parallel to the free surfaces). In both geometries the 
appropriate Green's function depends onR — R! because of homogeneity in those directions, mean- 
ing that only differentiation with respect to z and z' can give a non zero result for the "vector 
potentials" of Eq. (|l9|); in other words, only the component I — z of Af 171 is non zero. The equation 
to be solved for A 1 /" 1 is then: 

j^ m (z', z- R'-R) + ^rf m (z', z- R'-R) = e jr ,^f(2' R'-R)- (27) 

We have obtained the Green's function appropriate for these geometries, satisfying Eqs. (|l3|) and 
dl4|), using the method of Chishko (1989), who found the time dependent Green's function for 
a thin plate. The Green's functions are written as the infinite medium Green's function (G^) 
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minus a term (Hi m ) that satisfies the homogeneous equilibrium equation and that cancels the 
normal stresses produced by the infinite medium's Green's function at the free surfaces: 

G im (z,z';R-R')=G^l(x-x')-H im (z,z';R-R'). (28) 

This strategy differs somewhat from the one emplyed by Gosling and Willis (1994), who used a 
similar split not for the Green's function, but for a combination of gradients of the Green's function. 
The representation used for H im for the half space is the following: 

H im (z, z'; R-R') = - J dR"G { ^(z; R - R")a^) m {-z'- R" - R>) , (29) 

and a similar form for the thin plate: 

H im (z, z'; R-R') = ]T I dS'^iz', - i?'>i° )m (£ (a) - z; R" - R) , (30) 

a=l,2 

where the index (a) — 1,2 represents the free surfaces at z" — = ±h, and d5" = ±dR". We 
have introduced the surface Green's functions for the half space and the thin plate: they 

are the displacements in direction (i) produced by a localized impulse in direction (n) applied at 
the free surface. They satisfy, respectively, the homogeneous equilibrium equations: 

^a^^™^"^ = ° (31) 



Ci ^dx-dx- G ^ {z ^ n - n ' ] = ° (32) 



(£ = ±/i), and they are subject to the boundary conditions: 

a£ )n (z = 0;R-R') = -5 in S{R-R') (33) 
a^ n (^ b \C {a) ;R-R') = ±5 {am 6 m 5(R-R'). (34) 

It is possible to solve for the Fourier components (in the plane parallel to the free surface) of 
the Green's functions G^ {z; R - R') {g[*\z\k)), and g£J R - R!) (g$ 0,£|% 

g%\z\k) = JdR^G^{z-R) 



g%(z\k) = J dRe"-*G%{z;R), (35) 

and they are given in Appendices I and II. From Eqs. ( p9|) and (|30|), one obtains in Fourier space: 
4^ m (z',z\k) = -a^ n (z>\k)a^(-z\k) (36) 
a?^ m (z',z\k) = ^(±)aW n (z\m^ m (t-*\k), (37) 

with £ = ±h. Expressions for a^ h ^ m (z'\k) and <Jn} m (z — z'\k) are shown also in Appendix I, and 
for a^ p ^ m (z', £\k) in Appendix II. 
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Since 07° depends on (z' — z), only the term contributes to the determination of the 

"vector potentials" in Eq. fl27|). If a* m (z', z; R' — R) is the "vector potential" associated with cr^ m , 
namely 

da 
dx 



a% m {z', z- R'-R) = £jrs -i r (z', Z ; R! - R) , (38) 



then the "vector potential" A l s zm that satisfies Eq. ( p7f ) is given by: 

Al zm (z\z-R' -R) = -(£ + jL)a}™{z\z;R! -R) . (39) 

For the sake of brevity we will present the stress function tensor associated with of™ 1 , from which 
a™ and further ^4* zm can be determined by differentiation. The stress function tensor %™ is 
defined through a "double curl" : 

<J% m (z\ z; R'-R) = -e ipq e jln -^^{z' , 2; R' — R) , (40) 

a form that assures that f^ m , symmetric in (r/), is divergenceless in the indices i and j, as it 
should. It is easy to see that: 

a\ m {z', z; R'-R) = -e ipq -^(z', z; R' - R) . (41) 
In Fourier space the stress function tensor for the isotropic half space geometry is: 

— k(z-\-?') 

X z vv {z',z\k) = - ^ _ 1)fc4 i(7 2 - 2)[(7 2 + 1) + 2( 7 2 - l)kz]k n k u 

-il 2 ~ 1)[7 2 + (7 2 - Vkz + ( 7 2 + l)kz' + 2( 7 2 - l)k 2 zz']q v q u } 

^ e -k(z+z') 

X% v (z', Ak) = 2 1 2 { 1 2 - l)k 5 ^ 2 ^ 2 ~ l ^ q ^ q7 > kv + qvk ") 

+ ( 7 2 - 1)[1 - ( 7 2 - l)kz + ( 7 2 + l)kz' - 2( 7 2 - l)k 2 zz']kpq v q v 

-( 7 2 - 2)[( 7 2 + 1) - 2( 7 2 - l)kz]k k ri k w } , (42) 

where q v = e^gks, greek indices mean coordinates x or y, and 7 2 = cf/c 2 is the ratio between the 
longitudinal and transverse sound velocities in an isotropic medium. These formulae were obtained 
using the explicit forms for the Fourier components cr^ m (z', z\k) that are written in Appendix I, 
Eq. (|6l|), and solving for a stress function tensor that has only components x™u 7^ in order to 
reproduce <7^ m - Similarly for a plate: 

X ^(z',z\k) = E( ± ) 47 2 e (7 rr i ' )fc4 {(7 2 ~ !)[(! + (7 2 - l)k\Z - z\)(kz'fa - ktfa) 

-8{t - z)( n 2 + ( 7 2 - l)k\£ - Z\)(<f> 3 + k^ - kz'cp 2 )]q n q v 

+ ( 7 2 - 2)[a(£ - z)( 7 2 + ( 7 2 - ~ + (1 + (7 2 - - «|)&] W 
X ^(z',z|fc) = E( ± ) 47 2 (7 2 _ 1)fc5 i(7 2 - !)[(! ~ (7 2 ~ 1)^1^ ^ ^l)(03 + fcg0i ~ fcz> 2 ) 



-( 7 2 - 2)[ S (C ~ ^)(7 2 - (7 2 - 1)^1$ - *|)& + (1 - (7 2 - l)k\i - *\)<hV*f>h,K 
rj. s 9, 9 w cosh(fcz') sinh(fcz'), , , , , , . 

- *)tV - + )^(^ + ■ ( 43 ) 

5 Example: a screw dislocation perpendicular to the free 
surfaces of a half space and a thin plate 

Eshclby and Stroh (1951) obtained the elastic fields corresponding to a stationary screw dislocation 
perpendicular to the free surfaces of a thin plate. From their method it is easy to obtain the corre- 
sponding fields for a half space. As examples of the use of the general line integral representation 
for the stresses presented in this paper, Eq. (^21), we rederive these exact results. In the following 
we go through the main steps in the calculation. 

5.1 A screw dislocation in a half space 

Using the method in Eshelby and Stroh (1951), a simple calculation gives the displacement and 
stress fields, in polar coordinates, due to a screw dislocation lying along the z axis perpendicular 
to the free surface of a half space z > 0: 

U z (0) 
U e (R,z) 

a z e(R, z) 











-I 


00 dk 


2* Jo 


T e 


fib 




2ttR 




/ib 


\ dh 




'0 



-kz 



a re (R,z) = / dke- kz J 2 (kR), (44) 

27r Jo 

with J\(z) and J%{z) the Bessel functions of order one and two respectively, and /i the shear 
modulus. 

Now we describe how these results can be obtained as a special case of the general formula, Eq. 
(|22|). When applied to this case it reads: 

a Pq (x)=a$(x)+a${x), (45) 

where: 



a«(z,i?) = bC pqim e vl dz'a% m (z',z;-R) (46) 

Jo 

a$ (z, R) EE -bC pqmz [ dR' fj Af m {z' = 0, z; R' — R) , (47) 

J-c 

with — C a curve that closes the dislocation loop on the free surface, and £12 = —£21 = 1, en = 
£22 = 0. 
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First we focus on the first term a!pq (x) . For this term the segment on the free surface doesn't 
contribute due to the boundary condition: crf™(x' s , x) = 0. In Fourier space Eq. ( |46| ) becomes: 

POO 

o-W ( Z) -k) = bC pqim e vl / dz'a% n (z', z\k) . (48) 
Jo 

After doing the integration with the expressions in Appendix I, one obtains <Ji~z(z\k) = 0, and: 

a^(z\k) = -ib^. (49) 
Going back to real space, and in polar coordinates, we get 

cf^g(z\R) = /„fyf 2 / dk~~~T, — (sin sin a + cos cos a) , (50) 
(2t) J k 

where a is the polar angle associated with k, and 6 the angle of R with respect to some fixed axis. 
Defining <f> = a — 9, one gets effectively an integral over cj>, which is: 



<e\z\R) = J^J dk J dfe-*"""'™* = ^ ^ dkJ.ikR) = ^ , (51) 



i.e. it reproduces the first exact term of a z g(R,z) in Eq. (44). Similarly it can be shown that 
a zr (R, z) — 0. Also, this type of term is the only contributing to the calculation of a r e(R, z), and 
by a similar integration one gets the exact result in Eq. (44). 

Now we focus on the second term a^q(x). For this term, only the segment on the surface 
contributes. From expressions in Appendix I, we see that the "vector potential" A^ zz is zero at 



the free surface, and that Ag zS is given in Fourier space by: 



-kz 



Af 6 {z' = 0, z; k) = —kpe Sri k v . (52) 

Now, the integral § _ c dR'^Af 5 (x' , x) appearing in Eq. (^) is independent of the choice of curve 
— C on the free surface because the component z of the curl of A zzS with respect to the index j is 

zero on the free surface: indeed, ik^t-^pA^ °{z — 0, z\k) — 0. Thus, the expression for cr zS (x) in 
Eq. d47| ) becomes: 

afj(x) = - A f da J dR'- k J™ * e -^'-«l ^e Sn k v e~ kz , (53) 

with a the angle of k and <j> the angle between k and (R' — R). The curve C was chosen in a 
convenient way. First, it starts at R' = 0, which coincides with the position of the "vertical" 
segment of the dislocation. Then it continues in a semicircle around the point R (there \R' — R\ = 
R). After that, it continues on a straight line parallel to the direction of the vector R ( R! —R = XR, 
A G (1, oo)). In this way, the result from both of these segments adds up to the result: 



/>oo 

a^(R,z) = -^J dke- hz J x {kR), (54) 



which is the second term of the exact result in Eq. (J44J) , as it remained to be proven. 
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5.2 A screw dislocation in a plate 

Eshclby and Stroll (1951) obtained the exact displacement and stress fields, in polar coordinates, 
due to a screw dislocation along the z axis perpendicular to the free surfaces of a plate (the plate 
exists for —h < z < h): 



u z {9) 



b 



b f°° dk sinh(kz) „ , 
U e {R,z) = -— / -: T77TT J i( fci? ) 



2ir J k cosh(fc/i) 

in \ lib f°° cosh(fcz) 

<j ze (R,z) = - — — — — / dk — — — ■ Ji(fci?) 

2irR 2ir J Q cosn(fc/i) 

a re {R,z) = £ [°° dk^^-MkR) . (55) 
2n J cosh(kh) 



The verification that these results follow as a special case of Eq. (|2|) is analogous to the half-space 

\ zz 

l 



case: one needs the "vector potential" Ag zS evaluated at the free surfaces, which in Fourier space 



is, from Appendix I, 



A"S(+h ,-Pi- k P £5 » k » r | cosh(fcz) sinh(fcz) 

A [±h,z,k)~ 2~k^ {± cosh(kh) + S mh(kh) } - (56) 



The terms a^q (x) and a^j (x) from Eq. (|22|) become: 



h 



a$(z,R) = bC pqjm e VJ dz'af™(z', z; —R) (57) 



-ft 

pg v ' J - 2bL pqmZ{ ^ )2 J dk J^dR^e 2 fc2 CO sh(fc/ l ) • (58) 

The necessary integrations are done similarly as for the half space case (the curve C is the same) , 
and the only non zero components of the stress that are obtained are those given by Eq. (^5j). 
{ijpq (x) contributes the first term of a z e, and a^ q \x) the second term of (T z e). 



6 Discussion 

We have given a line integral representation for the stresses of an arbitrary dislocation loop in 
an arbitrary, not necessarily isotropic, three dimensional, finite elastic body. This representation 
is given in terms of a "vector potential" , for whose computation the Green's function for the 
elastic body in question must be known. This, in general, is not the case, but use of approximate 
expressions will yield approximations to the corresponding stresses. At any rate, we have provided 
explicit forms of the said "vector potentials" for a half space and for a thin plate. Previous results 
for the stresses due to a screw dislocation in a thin plate (Eshelby and Stroh, 1951), as well as 
easily obtained stresses for a screw dislocation perpendicular to the surface of a half space, were 
recovered as special cases. 

Gosling and Willis (1994) have provided an alternative line integral representation for the 
stress due to an arbitrary dislocation in an isotropic half-space. Our formulation is more general, 
admittedly at the price of introducing the Green's function for a finite elastic sample, a quantity 
that is in general hard to compute. 
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Appendix I: Green's function of a half space 

We write the Fourier components of the surface Green's function appropriate for a half space, 
g^j(z\k) (they arc the solutions in Fourier space of Eq. (|3l| ) subject to the boundary condition 
given by Eq. (||) : 

— kz 

;{7 2 + (7 2 -I)**} 
r {l + ( 7 2 -l)fcz} 

r{-l+(7 2 -l)^} 

:{2( 7 2 - l)q a q fi + [ 7 2 - ( 7 2 - l)kz]k a kp} , (59) 

where q a = e au k v is a vector perpendicular to k v (€12 = 1, £21 = — !)■ The Fourier components of 
crj° )m {z-z'\k) are: 



afXm = 


e -fcz 




2/i( 7 2 - 


i)k 


9 { ${z\k) = 


ikpe~ 


kz 


Mi 2 - 


l)k 




ikpe~ 


kz 


Mi 2 - 


l)k 




e -kz 


Mi 2 - 


l)k- 



o^{z-z'\k) = -^LLe-"'- '{ 7 2 + ( 7 2 -l)fc|z-z'|} 
af z z {z-z'\k) = -i^ e - fe l^'l{l + ( 7 2 -l)fc|^-/|} 



_ = _fl£ 7.|--' 4 

2 7 2 

f ")^. .'If. - ik P p -k\z- Z '\ , 1 ,.J 

2 7 2 fc 

a£) 5 (z-z'|fc) = ^ e -^-'l { l-( 7 2 _l)fc| z - z '|} 

ajf(z-z'|fc) = -f^Zjpe-kl*-. 'I { 7 2 g^ + [ 7 2 - ( 7 2 - l)fc|z - . 



(60) 



Doing the calculation in Eq. (p6[), one gets the Fourier components afj m (z', z\k): 
a^'(z',z\k) = ] { 7 2 + ( 7 2 - l)kz + ( 7 2 + l)fcz' + 2( 7 2 - l)fc 2 zz'} 



a^(z',z\k) = -^ e -^'){l + ( 7 2 -l)A : z-( 7 2 + l)A : z'-2( 7 2 -l)fc 2 zz'} 

-fc(z+z') 



a^(z',z\k) = ^ ' _ i)fc2 {( 7 2 - 2)[( 7 2 + 1) + 2( 7 2 - l)fcz]g a g/3 

+ (7 2 - 1)[(7 2 + 2) + 3( 7 2 - l)kz - ( 7 2 + l)kz' - 2( 7 2 - l)k 2 zz']k a kp} 
a£ s (z',z\k) = 0^e- k ^{l-(l 2 -l)kz+( 1 2 + l)kz'-2(^-l)ezz'} 
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e -k(z+z') 



af z {z\ z\k) = ' tfqgqp + [ 7 2 - ( 7 2 - 1)^ ~ (7 2 + l)kz' + 2( 7 2 - \)k 2 zz']k & kp} 



z\k) = zrm ™ ~ 2)[(7 + 1} ~ 2(7 " Wtaag/j + 

( 7 2 - 1)[2 7 2 + 1 - 3( 7 2 - l)fcz - ( 7 2 + l)fcz' + 2( 7 2 - l)k 2 zz']k s k a k p } 
+7 2 (7 2 ~ l)qs[k a qf3 + k q a ]} . (61) 

Finally, we derive the form of the Fourier components of the "vector potential" A% zz (z', z\k) and 
A¥ s (z' ', z|fc), necessary for the integration of Eq. Q). From Eqs. ([D]) and (^): 

/c ( £ | £ ' ) 

a£V, z|fc) = -ift,**, = 2jfc2 [ 7 2 + (7 2 - l)kz + ( 7 2 + l)kz' + 2( 7 2 - l)fc 2 zz'] , (62) 

and from Eq. (|39|): 

A^(z',z|fc) = -(^ + ±) a ^(z',z\k) = -^§l k z'e- k ^[l + ( 7 2 - l)kz] , (63) 



or A% zz {z' = 0, z|fc) = 0. Similarly: 



-k(z+z') 



af(z',z\k) = -iq v x S vl 3 = £ 2 2fc3 {l 2 qski3 

+ [1 - ( 7 2 - l)kz + ( 7 2 + l)kz' - 2( 7 2 - \)k 2 zz']k 5 qp} , (64) 

and from Eq. (|39|): 

Af\z\ z\k) = -(£ + ^KV, z|fc) = '""IP } rf k W + 2 ^ ~ & - ^)kz]kz'k sqp } , (65) 
or A zz& {z' = 0,z\k) = exp(-kz)k fj q s /k 2 . 

Appendix II: Green's function of a thin plate 

We write the Fourier components of the surface Green's function appropriate for a thin plate, 
g\f{z\k) (they are the solution to Eqs. (ph and (63)): 



s « = M^—Tjk {l2Mzlc) ~ { ^ - i)Ms(*IO + (t 2 -i)KM*\t)} 



9$(*,t\k) = 1)fc2 {(7 2 - l)kzM*\Z) (7 2 - l)fcto(*|fl - 04(^|C)> 

= - 4M7 f^ 1)fc2 {(7 2 - i)fc^i(z|0 - (7 2 - l)fcz02(z|O - &(s|0} 

^S^'fl*) = ^p{(7 2 - 1)^4(4) - (7 2 - 1)^3(4) + 7 2 *(^)} 
(+] q a qi3 f cosh(fcz) sinh(fcz) 

1 j 2^ 3 sinh(fc£) cosh(fc£) J ' 1 ' 
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where £ = ±/i, and the functions </>i(z|£) to </>4(z|£) are defined in the following way: 
01 ( z \0 = t~ cosh(fcz) cosh(fc£) + — sinh(fcz) sinh(&£) 
4>2{ z \€) = T" sinh(fcz) sinh(fc^) + — cosh(kz) cosh(fc£) 

Os a 

03(z|£) = T" cosh(fcz) sinh(fc£) + — sinh(fcz) cosh(fc£) 

04(z|£) = — sinh(fcz) cosh(fc£) + cosh(fcz) sinh(fc£) , (67) 
Os o a 

with: 

(5 S = smh(kh) cosh(fcft-) + fc/i 
5 a = smh(kh) cosh(fc/i) — kh . (68) 

Doing the calculation in Eq. (|37[), one gets the Fourier components cr^ m (z', z|fc): 
af/(z',z|fc) = ^(±)l^^{(l + ( 7 2 ^l)A ; |^z|)(fcz> 4 (z'|a"fc^3) 

- ^)(7 2 + (7 2 - ~ z|)(& + - fc*'&)} 
of/(z',z|fc) = -^(±)!^^ {s (e-z)( 7 2 + ( 7 2 -l)A ; |e-z|)(A : z'03(^IO-fc^4) 

+ (1 + ( 7 2 - l)fc|£ - z|)(0 4 + fcz'fc - k^ 2 )} 
a^(z',z\k) = -EW 47 ^r' 1) {^(7 2 -2)x 

" ^)(7 2 + (7 2 " " *|)& + (1 + (7 2 " " *|)^i] 
+fc Q A: /3 ( 7 2 - l)[s(£ - z)( 7 2 + ( 7 2 - l)k\£ - z|)(0 3 + ~ 
+ (1 + ( 7 2 - l)fc|£ - z|)(20x + - 

a? z S (z',z\k) = ]>>) lhe ~^ k A {(1 ~ (7 2 ~ l)*|g - z\)(<h + fcgfr - kz'fo) 



-*(£ - z){kz'<j>, - fc£0 3 )( 7 2 - ( 7 2 - l)fc|C - z\)} 

<J^(z',z\k) = ^(±)f-^i{fc 5fc/3 [(l-( 7 2 -l)fc|e-z|)(fc^4-fc^3) 



_ s (£ _ 2) ^ 4 + fcz ' 01 _ fc^ 2 )( 7 2 - (7 2 - mt ~ z \)} 

2 sinh(fcz') cosh(fcz') 

-(±)ff, w .(f - *) 7 [— ^ + ^TT^y]} 

g-*|£-«| 

a^(z',z\k) = E(±) 472(/y2 _ 1)fc3 {^ fc « fc M7 2 - 1) x 

[(1 - ( 7 2 - l)fc|£ - Z|)(0 3 + fcz'0 2 - fc^l) + 

«(£ - z )(7 2 - (7 2 - - *|)(2& + fcz>4 - fc^ 3 )] + ik S q a qp(l 2 - 2) x 
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[(1 - ( 7 2 - - z\)h + *(£ - z)(7 2 - (7 2 " " z\)<h] 
.fa-V /t \ 2/ 2 lU cosh(fcz / ) sinh(fcz') 



(69) 



with s(£ — z) = sign(£ — z). Finally, we write the Fourier components of the "vector potential" 
evaluated at the free surfaces, necessary for the integration of Eq. (J5q) : 



Af z {z' = ±h,z\k) = 

A" s ( s /-+h X \P\ - k p qs coshjkz) sinh(fcz) 
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